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Macroscopic quantum tunneling is described using the master equation for the reduced Wigner 
function of an open quantum system at zero temperature. Our model consists of a particle trapped 
\^ • in a cubic potential interacting with an environment characterized by dissipative and normal and 

f"*) | anomalous diffusion coefficients. A representation based on the energy eigenfunctions of the isolated 

system, i. e. the system uncoupled to the environment, is used to write the reduced Wigner func- 
£N| ■ tion, and the master equation becomes simpler in that representation. The energy eigenfunctions 

computed in a WKB approximation incorporate the tunneling effect of the isolated system and the 
effect of the environment is described by an equation that it is in many ways similar to a Fokker- 
Planck equation. Decoherence is easily identified from the master equation and we find that when 
the decoherence time is much shorter than the tunneling time the master equation can be approxi- 
mated by a Kramers like equation describing thermal activation due to the zero point fluctuations 
of the quantum environment. The effect of anomalous diffusion can be dealt with perturbatively 
and its overall effect is to inhibit tunneling. 

> 

t""— , PACS numbers: 03.65.Yz, 03.65.Sq, 05.40.Ca 

o . 

. I. INTRODUCTION 

o 

The observation of quantum tunneling effects is now possible in some macroscopic quantum variables such as the 
flux quantum transitions in a superconducting quantum interference device, or the decay of a zero- voltage state in 
a current-biased Josephson junction [ij, |2j, |3|, |^| . Macroscopic quantum systems pertains to the boundary between 
quantum and classical physics. These systems are modelled by open quantum systems which are characterized by a 
distinguished subsystem, the "system" for short, interacting with an environment. The environment acts as a source 
of dissipation and noise for the system and produces quantum decoherence which generally inhibits tunneling 0, |(| . 
3 There is a vast literature on this subject, see Refs. 0,H,Ej for comprehensive reviews. 

Most work on macroscopic quantum tunneling is based on imaginary time formalisms such as the Euclidean func- 
tional techniques which have been introduced in the classical field of noise-activated escape from a metaestable state 
QU, or the instanton approach introduced for quantum mechanical tunneling or for vacuum decay in field theory 
|l lL Il2l IT3I ItH IT5L ITfjl ] . These techniques are specially suited for equilibrium or near equilibrium situations, but are 
£3 . difficult to generalize to non equilibrium situations. 

There are theoretical and practical reasons for a formalism of non equilibrium macroscopic quantum tunneling. On 
the theoretical side dissipation as well as decoherence are only truly understood in a dynamical real time formalism. 
In the classical context thermal activation from metaestable states is well understood since Kramers ^3 m terms 
of the dynamical Fokker-Planck transport equation, where the roles of dissipation and noise and their inter-relations 
are known. On the other hand, an open quantum system may be described by a dynamical equation for the reduced 
density matrix, the so-called master equation, or the equivalent equation for the reduced Wigner function which has 
many similarities to the Fokker-Planck equation. However, at present no compelling derivation of the tunneling rate 
is available in this dynamical framework, that might be compared to the instanton approach for equilibrium systems. 
Consequently, the effect of dissipation, noise and decoherence on tunneling and their inter-connections is not yet 
fully understood. On the practical side out of equilibrium macroscopic quantum tunneling is becoming necessary 
to understand arrays of Josephson junctions, or time-dependent traps for cold atoms which are proposed for storing 
quantum information in future quantum computers [la . Il9t l2fl l2l| , or to understand first order phase transitions in 
cosmology |23.l23|. 

In recent years we have considered different scenarios in which metaestable quantum systems are described by the 
master equation for the reduced Wiener function. By using techniques similar to those used for thermal activation 
processes on metastable states |l7l |24| it was possible to compute the effect of the environment on the quantum decay 
probability. This was used in some semiclassical cosmological scenarios for noise induced inflation 25] due to the back 
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reaction of the inflaton field, in the context of stochastic semiclassical gravity prl I27I l28l I29I l30 | : see Refs. [3ll l3^| 
for reviews on this subject. It was also used for bubble nucleation in quantum field theory, where the system was 
described by the homo gen eous mode of the field of bubble size and the environment was played by the inhomogeneous 
modes of the field |33l l34j . and on some simple open quantum systems coupled linearly to a continuum of harmonic 
oscillators at zero temperature |35j| . But in all these problems only the contribution to tunneling due to activation 
was considered. One of the purposes of this paper is to clarify when this approximation is valid. 

In this paper we propose a formulation of macroscopic quantum tunneling using the master equation for the reduced 
Wigner function in which both the pure quantum tunneling effect and the environment are taken into account. This 
is possible by the introduction of a representation of the reduced Wigner function based on the energy eigenfunctions 
of the isolated system, i. e. the system not coupled to the environment. This representation is useful in a way 
somewhat analogous to the way the energy representation is useful in the Schrodinger equation. The key to this 
result is that quantum tunneling is already encoded in the energy eigenfunctions, which we can compute in a WKB 
approximation. It is quite remarkable that in this representation the master equation can be solved analytically under 
certain approximations. 

In order to have a working model in a form as simple as possible, but that captures the main physics of the problem, 
we use constant dissipation and normal and anomalous diffusion coefficients to describe the effect of the environment. 
These coefficients can be deduced from microscopic physics: they take constant values when the environment is made 
by an Ohmic distribution of harmonic oscillators weakly coupled in thermal equilibrium at high temperature; but at 
zero temperature they are time dependent |35l l3(i l37| . Thus the model studied here may be seen as a toy model at 
low temperature, generally valid at long time scales only. 

Master equations play also an important role in elucidating the emergence of classicality in open quantum systems 
as a result of their interaction with an environment. In fact, as the master equation gives the quantum evolution 
of initial states, defined by the reduced Wigner function at some initial time, it has been of great help to study 
decoherence. In particular, it has been used to clarify the way in which the environment selects a small set of 
states of the system which are relatively stable by this interaction, the so-called pointer states, whereas the coherent 
superposition of the remaining states are rapidly destroyed by decoherence |38t 13^. I40L l4ll li^ . Using large scale 
numerical simulations the effect of the interaction with the environment on coherent tunneling has been analyzed 
in the framework of an open quantum system that is classically chaotic: a harmonically driven quartic double well 
[43l lil ] . More recently 0] tunneling in a simple double well potential has been numerically simulated using the 
master equation at high temperature as well as at zero temperature. It is found that at zero temperature tunneling is 
inhibited by the environment that produces decoherence nevertheless at large time scales tunneling is still possible by 
an activation-like process due to the zero point fluctuations of the quantum environment. In the model we consider in 
this paper, which involves a cubic potential, we find a very similar behavior when the decoherence time is very short. 

This paper is organized as follows. In the next two sections ^ and IIIII we review the theory of tunneling in 
closed systems and introduce the energy representation for Wigner functions. This extended review is necessary both 
to establish our conventions and to recall specific results which are central to the main argument. In section llVl we 
introduce the environment and write the master equation for the reduced Wigner function of the open quantum system, 
discuss the different dissipation and diffusion coefficients and derive the master equation in the energy representation. 
In section this master equation is explicitly computed and the different time scales of the problem are discussed. 
In section IVT1 under the assumption of strong decoherence we analytically solve the master equation, which is reduced 
to a Kramers-like equation, and then perturbatively compute the effect of anomalous diffusion on tunneling. Finally, 
in Section IVIII we briefly summarize our results. In the Appendices we provide additional technical details. 

II. TUNNELING IN QUANTUM MECHANICS 

In this section we review the WKB method to tunneling in quantum mechanics. The energy eigenfunctions in the 
WKB approximation we obtain will play an important role in the energy representation of the Wigner function that 
will be introduced latter. 



A. The system 

We begin with the simple closed quantum mechanical system formed by a particle of mass M in one dimension 
described by a Hamiltonian 



(1) 
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FIG. 1: A schematic plot of the potential. For an energy E below the barrier there will be three classical turning points, also 
shown. 

with a potential U given by 

U(x) = hm 2 x 2 ~^x 3 , (2) 

for small values of the coordinate x. This is a fairly general potential for a tunneling system, it is the basic clement 
in the dashboard potential, which is a very good model for a flux trapped in a superconducting quantum interference 
device (SQUID), or a single Josephson junction biased by a fixed external current 0,0,00] For technical reasons, 
it is convenient to assume that for large x the potential flattens out and takes the value U (x) — —Urx, both negative 
and constant. The tunneling process ought to be independent of the form of the potential this far away from the 
potential barrier. We present a sketch of this potential in Fig. 1. 

There is one classically stable point at x = 0, and one unstable point x — x s — 2MCIq/X, corresponding to an 
energy e s = 2M 3 ft[j/(3A 2 ). The curvature of the potential is U"(0) = MQ 2 , at x = 0, and -U"(0) at x s . The other 
point at which U (x) = IS X — X ' Qxit 

= (3/2)x s . For x 3> x ex u the potential flattens out and is constant. 

B. The WKB approximation 

If we assume that the particle is trapped in the potential well, that is in its false ground state or false vacuum, 
the tunneling probability can be computed in this simple problem in many ways. One of the most efficient is the 
instanton method which reduces to the computation of the "bounce solution" . The most attractive aspect of this 
computation is that it can be easily extended to field theory where the tunneling probability is then interpreted as 
the probability per unit time and volume to nucleate a bubble of the true vacuum phase. The rate for quantum 
tunneling is T c i osec i = (f2o/27r)a 9 cxp(— Sb/K)i where Sb is the action for the "bounce" (or instanton), namely the 
solution to the classical equations of motion which interpolates between x = and x = x ex it in imaginary time 
Sb = 2 Jq"^ dx ^J2MU (x) : and the prefactor a q = (1207rS's/?i) 1 ^ 2 . Our expression for the potential is so simple 
that the above integral can be computed explicitly: S's/ft = 18e s /(5eo)i where £o = \K^o is the zero point energy of 
a harmonic oscillator with frequency Qq. 

Here, however, we will concentrate on a real time approach by expanding the false vacuum state as a linear 
combination of true eigenstates of the Hamiltonian. To the required accuracy, it is enough to work with the WKB 
approximations to the true eigenfunctions; see for instance Refs. 0>E3- The instanton method reviewed in the 
previous paragraph can, in fact, be easily justified by this semiclassical approximation. Here we explain in some 
detail this standard procedure to obtain the eigenfunctions by matching the WKB solutions in the different regions 
of the potential. These solutions will play a crucial role in the energy representation for the Wigner functions to be 
introduced latter. 
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Let < E < e s be the energy of the particle in the potential well, and ipE the corresponding eigenfunction. The 
Schrodinger equation is 

~' r '' 2 ^E + U(x)^ E =E^ E . (3) 



2M dx 2 
Let us define 

p(x) = ^/2M\U(x)-E\, (4) 
and the integral S(x, y) (note the order in the integration limits) 

S(x,y)= f dx'p(x'). (5) 
Jy 

The WKB solutions are obtained from these elements. We have to match the WKB solutions in the different 
regions across the potential function. The details of this calculation are given in Appendix A. The WKB solution 
i/je for energies in the range < E < e s is given by Eq. I|A14|) . where xl < xr < x out are the three classical 
turning points for the cubic potential J2J); see Fig. 1. The normalization constant Ke in Eq. (|A14|I is obtained by 
imposing the continuous normalization of the cigcnfunctions given in Eq. llATojl and it is given in Eq. (|A23|) . Of 
particular relevance is the value of the eigenfunction tjjE(x) at values x 3> x ou t- This gives the main contribution to 
the continuous normalization integral. The value of the eigenfunction at i > x ou t, as computed in Appendix A, is 



where the phase Se is introduced in Eqs. (|A24I) and Poa{E) is defined by Eq. (0} when x ^> x out ; see also Eq. <|A17|) . 

We are interested in the details of the eigenfunctions near the false vacuum state, since we will be dealing with 
tunneling from vacuum. Thus, in the remaining of this section we give explicitly the values of the normalization con- 
stant Ke and the phase shifts 8e near this vacuum state. Therefore let us impose the Bohr-Sommerfeld quantization 
condition (|A45|) and let Eq be the corresponding lowest energy, that is, n = in Eq. As we will see in the 

next subsection this defines the false vacuum energy. Expanding the integral in Eq. Q around Eg we find that close 
to the lowest energy value 



S(x R ,x L )~?f-T(E-E ), (7) 



where r is defined by 




(8) 



U (x) - E ' 



Thus cos 2 (S (xr,xl) /K) ~ (t 2 /h 2 ) (E — Eq) 2 , and evaluating the right hand side of (|A26|) at Eq, we conclude that 
K E has poles at the complex energies 

E± = E ± is, e = A e -2So(x „ t) tf R )/fi / 9 ) 
4r 

which is in agreement with the standard result 48]. To simplify the notation let us call Sq = So (x out , xr) and 
/o = f (E Q ) + 7r/4, then we have from Eqs. HA21jl and l|A22j) that the functions A(E) and B(E) for E near E 
are: A (E) = (t/H) exp(S /h)[F_(E) + F + (E)} and B (E) = (-ir/h) exp(S /h)[F^(E) - F+(E)}, where F-(E) = 
exp(if a /h) (E — and F + (E) = exp(—if /H) (E — E + ). Notice that neither A nor B vanish at E±. Finally from 
Eq. (|A23p we can write the normalization constant near the false vacuum energy, as 

2 M e _ AMe 2 e 2Sa ' h 

E ~^(E-E f + e 2 ~ ttH 2 (E-EJ){E-E + Y ( ] 

and from Eqs. 1A24I) the phase shifts are 



£ 2 e 2S /h if , h E-E- 



V (E-E ) 2 + e 2 K ' VE-E+ { ' 
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Equations ©, i|10|) and are the main results of this section. We notice, in particular, the poles of the norm 
and the phase shifts at E± near the false vacuum energy. The strong dependence on the energy of these functions 
near the false ground energy will play an important role in the next sections. We will need, in particular, the phase 
shifts derivatives which are given by: 



d5 E _ -i f 1 1 
~dE ~ T V E-E- ~ E-i 



(12) 



C. The false vacuum 



Before we start with the computation of the tunneling rate we have to define what we mean by the decaying state, 
all the wave functions we considered so far are true stationary states and, obviously, show no decay whatsoever. We 
need to confine initially the particle into the potential well in its lowest energy. To this end, we introduce an auxiliary 
potential U aux which agrees with U up to x s (where the true potential reaches its maximum value) and increases 
thereafter. We may assume that the growth of £/ aux is as fast as necessary to justify the approximations below; the 
tunneling rate is insensitive to the details of U aux beyond x s . Thus, we define the decaying state ipo as the ground 
state of a particle confined by ?7 a ux |49| . 

It is obvious from the form of the WKB solutions that ipo agrees with ipE up to x s , i. e. ipo(x) = ipE (x) for 
x < x s , where Eq is the Bohr-Sommerfeld ground state energy for the auxiliary potential C/ aux , which corresponds to 
n = in the condition l|A15|) . Beyond x s , ipo will decay rapidly to zero, unlike ipE - Like any other wave function, ipo 
admits a development in the complete base of energy eigenfunctions ipE, as 

V>o (x) = J dE C e ^e (x) , (13) 

where due to our normalization the Fourier coefficients are given by 

C E = J dx ip E (x) ipo (x) . (14) 

To find these coefficients, we observe that ipo (x) is a solution to the Schrodinger equation with the auxiliary potential 

-h 2 d 2 

ipo + E4ux (x) ip = E ipo. (15) 



2M dx 2 



Let us add to both sides of this equation the term [U (x) — J7 a ux (^)]^o and then multiply both sides by tpE (x) and 
integrate to obtain 

/•OO 

(E-E )C E = - dxip E (x)[U am (x)-U(x)}iP (x). (16) 

J x s 

An important consideration is that ipo (x) is a smooth function (as opposed to a distribution), and, unlike tpE it is 
normalizable, so Ce must also be smooth. This means that it is allowable to assume E ^ Eq in Eq. (|16|l : Ce can 
then be found by analytical continuation. To estimate the right hand side of Eq. (|16(l . let us introduce; cf. Eq. (J2J, 



Paux {x ) = ^2M\U am (x)-E \. (17) 
To the right of x s we may use the WKB approximation with the decaying solution into the forbidden region to write 

ipo(x) = ^ (x s )exp (~^J Paux(y)dy^j. (18) 

On the other hand, ip E (x) is given by Eq. (|A12|I in Appendix A. If E is close to E , then Eq. (7J applies, and we 
may write 



ip E (x) - 2K E 



(E - E ) F + (x s ,x R ) exp Q J p (y) dy^j + (x s , x R ) exp (-^ J 



P (y) dy 



(19) 
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Substituting the two previous expressions into the right hand side of Eq. (|16fl we see that we have to compute the 
two following integrals, 

J±= f da [U mx {x)-U(x)]exp(~ f [p aux (y) ± p (y)} dy ) . (20) 



h 

The integral, J_, is dominated by the region near the lower limit, where p aux (x) is close to p (x) and we can write 



from where we obtain 



M 



-{E-E )[dx exp(-i£ 



J-=H ^4r^-(E-E ) I dx exp(-i I \p awc (y)-p(y)]dy) ) (21) 



where the remaining integral is made negligible by an appropriate choice of t/ aU x- For the other integral, J+, we see 
that the corresponding exponential factor in Eq. (1201) decays faster than the exponential factor of J_, so that the 
region which effectively contributes to the integral is narrower. Since the pre-exponential factor vanishes at the lower 
limit, we find J+ ~ 0. Finally, putting all these pieces together into the right hand side of Eq. (|16[) we get to leading 
order, 



(E - E ) 

whose solution, assumed smooth, is 



OTJ ( r ) 

C E + 2K E ^ (x s ) T\j -j^ ll F + (x s ,x R ) 



0. 



OTJ (r ) 

C E = -2K E ^ ( Xs )t^^^-F + (x s ,x R ) . (22) 

We note that Ce is independent of the choice of U aux beyond x s , as it should. 

Thus, we have found the false vacuum wave function in terms of the energy eigenfunctions of the original problem. 
The false ground state is a superposition of energy eigenstates which are fine tuned in such a way as to produce 
destructive interference outside the potential well. Notice that C E , because of the factor Ke in Eq. I|22|l . peaks near 
the energy of the false ground state, and has a strong dependence on the energy near this ground state energy. A 
good approximation for C E is given by 

C| = n (E- E + )(E - EJ) ' (23) 

D. Tunneling from the false vacuum 

Let us now compute the tunneling rate assuming that the particle is described initially by the false ground state 
-00- At time t, we have 

*/>(x,t) = J ' dEe- iEt ' h C E ^E{x), (24) 

The persistence amplitude is 

p (t) = / dx V>o (x) ip {x, t)= I dE e-' lEt/h C 2 E . (25) 



With the value of C E given by Eq. I|23|) p(0) = 1. To perform the integration we can close the contour of integration 

^2 



in the complex E plane adding an arc at infinity, whereby we pick up the pole _E„ in C' E , therefore p(t) goes like 



p(t) ~ exp 



^exp \ —\ S {x ou t,x R/ 



(26) 



(no prefactor) provided t is not too large. The tunneling rate for this closed system, T c i osec i, may be defined from the 
persistence probability p 2 {t) ~ exp(— T c i ose dt), so that T c i ose d — (l/2r) exp(— 2Sq/K), which agrees with the result of 
the bounce solution. Note that if we take the classical lowest energy E = 0, then x R — xl — 0, x ou t — x ex u, and 
Sb = 2S(x eX it,0), but here is the action corresponding to a particle with false vacuum energy Eq, which differs 
from zero, consequently it differs from Sb/2. This difference is accounted for by the prefactor a q in the instanton 
result. An equivalent way of deriving this result is to estimate the integral by a stationary phase approximation. 
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III. WIGNER FUNCTION AND ENERGY REPRESENTATION 



An alternative description of a quantum s yste m is given by the Wigner function in phase space, which is defined by 
an integral transform of the density matrix |5fl l51j . The Wigner function for a system described by a wave function 
Tp{x) is 



(27) 



where the sign convention is chosen so that a momentum eigenstate tp Po (x) 
(l/2irh)5 (p — po). Moreover, it satisfies 



x / h /y/2wh becomes W Po (x,p) 



dpW(x,p) = \ip (x)\ 2 , 



dx W (x,p) = 



e -ipx/H 

dx — ip (x) 



(28) 



and it is normalized so that J J dx dp W(x,p) = 1. Thus the Wigner function is similar in some ways to a distribution 
function in phase space, it is real but, unlike a true distribution function, it is not positive defined; this is a feature 
connected to the quantum nature of the system it describes. 
The Schrodinger equation for the wave function tp 1 



2M^ + U{x) ^ = th di^ 



(29) 



translates into a dynamical equation for the Wigner function, which is easily derived. In fact, by taking the time 
derivative of H27J) . using the Schrodinger equation (|29|l . and integrating by parts we have 



0_ 

dt 



W(x,p) 



*(-S)K-S)-"("!)k('+!)}- 



f dy Apy/n { 


f-ihp\ 




1 2ttH \ 


V M ) 


dx - 



For the cubic potential © we have U (x - y/2) - U (x + y/2) = -M^xy + (\/2)x 2 y + (A/24)y 3 and, noting that 



lie 



ivy in _ —ifrQ pe wy/ h a nd y 3 e lpy ^ h = ih 3 d p e lpy ' h , we get the equation for the Wigner function 



§-W(x, P ) 



n't \ 9 P d \ 2 d 3 
{X> dp Mdx + 24 dp 3 



W(x,p), 



(30) 



which may be interpreted as a quantum transport equation. The first two terms on the right hand side are just the 
classical Liouville terms for a distribution function, the term with the three momentum derivatives is responsible for 
the quantum tunneling behavior of the Wigner function in our problem. A theorem by Pawula ,52!] states that a 
transport equation should have up to second order derivatives at most, or else an infinite Kramers-Moyal expansion, 
for non-negative solutions W(x, p, t) to exist. The above equation for the Wigner function circumvents the implications 
of the theorem since it need not be everywhere-positive. Even if we have an everywhere-positive Gaussian Wigner 
function at the initial time, the evolution generated by an equation such as Eq. (|r?()|> will not keep it everywhere- 
positive. Thus, here we see the essential role played by the non-positivity of the Wigner function in a genuinely 
quantum aspect such as tunneling. 



A. The energy representation 

Given that a wave function ip can be represented in terms of the energy eigenfunctions ipE, defined by Eq. Q, as 

^ (x) = [ dE Ce^e (x) , (31) 



we can introduce a corresponding representation for W{x,p) in terms of a base of functions We 1 e 2 { x >p) m phase 
space defined by 



(32) 
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Then W (x,p) can be written as 

W(x,p) = J dE x dE 2 C ElE2 W ElE2 (x,p) , (33) 
where, in this case, we have Ce 1 e 2 — Ce^C e . On the other hand from the definition of W ElE2 (x,p) we can write 



lE2 (x,p) W E , E , (x,p) = / ^§ {ip El (x - I) xpE 2 (x + |) fe; (x - |) i) B > 2 (x+ I)} 



where the p integration has been performed. If we now call z — x — y/2, z' = x + y/2; then dxdy — dzdz', and 

/ ^W ElE Jx,p)W ElE ,(x,p) = ^5(E 1 -E[)5(E 2 -E! 2 ), (34) 

which gives the orthogonality properties of the functions W ElE2 . This suggests that any Wigner function may be 
written in this basis as 

W (x,p, t) = J dE x dE 2 C ElE2 (t) W ElE2 (x,p) . (35) 

We call this the energy representation of the Wigner function. In this representation, the master equation or the 
quantum transport equation l|30|) is very simple 

of ElE2 (t) = ^~{Ei- E 2 ) C ElE2 (t) , (36) 

as one can easily verify. One can give an alternative derivation of the tunneling rate from this equation, by taking 
the initial condition for the Wigner function which corresponds to the false vacuum. 

B. Tunneling in the energy representation: closed system 

Let us compute here the tunneling rate from the false vacuum for our closed quantum system. We assume that 
our particle at t = is trapped into the well of the potential J3J) in the false ground state with the energy Eq, i. e. 
the ground state of the auxiliary potential U aux introduced in Section III CI We know from that section that the wave 
function i/jq of this state can be expressed in terms of the eigenfunctions i/je by Eq. (|13fl with the coefficients C E 
given by Eq. I|22|) . or Eq. (|23|) . In terms of the reduced Wigner function, which we may call Wo(x,p), this state is 
easily described in the energy representation by the coefficients C ElE2 (Q) — C El (0)C E2 (0), where C E (0) is given 
by Eq. (|23(l . Because the dynamics of the quantum transport equation is trivial in the energy representation (|36fl the 
time dependence of the coefficients C ElE2 (t) is simply 

C ElE2 (t) = e-i^- E ^C ElE2 (0) . (37) 
Thus, according to Eq. I|35|l . the Wigner function at any time is 

W(x,p,t) = J dEidE 2 e~^ El ~ E2)t C El (0) Cb 2 (0)W ElE2 (x,p) . (38) 

From this we can compute, in particular, the probability of finding the particle at the false vacuum at any time. 
In terms of the false vacuum Wigner function and the Wigner function of the tunneling system we may define that 
probability as 

p 2 (t) = 2irh J dxdpW (x,p)W(x,p,t). (39) 

This equation can be used in a closed as well as in an open system. For the closed system of Section HP where the 
state is described by the wave function ip of Eq. 124JI and the false vacuum is described by the wave function i/jq 
of Eq. (|13|) . the square of the persistence amplitude 1251) is given, in fact, by Eq. I)39|l when the definition of the 
Wigner function, i. e. Eq. I|27l) . is used. For the open system the quantum state is not described by a pure state 
and, in general, the Wigner function W(x,p, t) can be written as W = ^2 i PiWi where pi is the probability of finding 
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the system in the state <j>i and W$ is the Wigner function for the state fa. The definition l|39|l leads in this case to 
P 2 {t) = YliPi I ("00 1 <^>i) 1 2 ^ which is indeed the probability of finding the system in the state i/V 
When the energy representation is used Eq. iffify becomes 

2 (t) = ( dE x dE 2 e-^-^Ci, (0) Cl (0) . (40) 



P 

To compute p 2 (t) we shall use the stationary phase approximation. The idea is that the integration paths for E\ and 
E2 may be deformed simultaneously in such a way that the integrand comes to be dominated by Gaussian peaks. For 
late times it is enough to seek the stationary points of \{E\ — E2). In principle, we could include C% and C% as fast 
varying components of the integrand, but these functions are really fast varying in the vicinity of E- and E + . Thus, 
the stationary phase condition leads to p 2 ~ exp(— 2te/H) in agreement with the persistence amplitude of Eq. J2HJ. 
The tunneling rate is T c i ose d = 2e/H. 

It is often convenient to give the tunneling rates in terms of an equivalent thermal activation problem under a 
potential barrier of hight e s with a certain effective escape temperature T esc . This escape temperature is defined Q 
from a given tunneling rate, T, by the equation 

r = _L e -s,/k B T eac _ ( 41 ) 

For the closed system, either the WKB approximation or the instanton method yield the same T esc with a barrier 
penetrability A, defined by V closed — (1/2t) exp(— A) (which gives A = (2/H)Sq in the WKB approximation); see Eq. 
(|D7jl in Appendix D. 

In the following sections we will use the energy representation of the Wigner function to compute the tunneling 
rate in a more complex problem involving coupling to an environment. The dynamics of the transport equation in the 
energy representation is simpler than in the phase space representation and the initial condition is given in terms of 
the coefficients (1221) which we have already computed. The task would be more difficult starting from the transport 
equation in phase space, such as Eq. (|3L)[) . since the third derivative term makes the solution of the equation very 
complicated. One has to resort to methods such as those based on matrix continued fractions in order to compute 
decay rates from master equations for open quantum systems with third order derivative terms |52l l53l l55| . 



IV. THE OPEN QUANTUM SYSTEM 

Now we assume that our system of interest is coupled to an environment. As emphasized by Caldeira and Leggett 
any quantum macroscopic system can be modelled by an open quantum system by adjusting the coupling of the 
system and environment variables and by choosing appropriate potentials. One of the main effects of the environment 
is to induce decoherence into the system which is a basic ingredient into the quantum to classical transition 0, l38l 

mm mm. 

The standard way in which the environment is introduced is to assume that the system is weakly coupled to a 
continuum set of harmonic oscillators, with a certain frequency distribution. These oscillators represent degrees of 
freedom to which some suitable variables of the quantum system are coupled. One usually further assumes that the 
environment is in thermal equilibrium and that the whole system-environment is described by the direct product 
of the density matrices of the system and the environment at the initial time, so that there are no initial system- 
environment correlations. The macroscopic quantum system is then described by the reduced density matrix, or 
equivalently, by the reduced Wigner function of the open quantum system. This latter function is defined from the 
system-environment Wigner function after integration of the environment variables. 

In order to have a working model in a form as simple as possible, but that captures the main effect of the environment, 
we will assume that the reduced Wigner function, which we still call W(x,p), satisfies the following dynamical equation, 



§-W(x,p) 



U » 7JZ " T7^Z + + «Z IP + lMa 2 - + HA- 



dp M dx 24 dp 3 dp \ dp dx 



W(x,p), (42) 



where 7 which has units of inverse time is the dissipation coefficient, and a 2 and A are, respectively, the normal and 
anomalous diffusion coefficients. The last three terms of this equation represent the effect of the environment: the first 
describes the dissipation produced into the system and the other two are the diffusion or noise terms. An interesting 
limit, the so-called weak dissipation limit, is obtained when 7 — > 0, so that there is no dissipation, but the coefficient 
7<7 2 is kept fixed. We will generally refer to equation (|42|l as the quantum Kramers equation, or alternatively, as the 
quantum transport equation. This equation reduces to a classical Fokker-Planck transport equation when h = 0: it 
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becomes Kramer's equation |24| for a statistical system coupled to a thermal bath and has the right stationary 
solutions. 

This equation can be derived 

laE EE El El El Ej 

assuming the so-called Ohmic distribution for the frequencies 
of the harmonic oscillators. In the high temperature limit, 7 is constant, a 2 — ksT, and A cx (fc^T)" 1 , where ks is 
Boltzmann's constant and T the bath temperature. Thus at high temperature one can generally ignore the anomalous 
diffusion term. In the low temperature limit, however, the master equation for the reduced Wigner function involves 
time dependent dissipation and diffusion coefficients. Typically the dissipation coefficient "/(t) starts with a zero 
value and after a short transient time, after which the system and the environment become correlated, it reaches a 
constant value; the normal diffusion coefficient a 2 (t) starts also with a zero value it reaches a maximum and after the 
short transient time it undergoes a mild oscillatory behavior until at time scales t 3> reaches a constant positive 
asymptotic value; the anomalous diffusion coefficient has a similar qualitative behavior but its asymptotic large time 
value is negative and depends on the cut off frequency. To be specific EE EE El EH ; at large time scales the normal 
diffusion coefficient becomes a 2 ~ ^M7o, and the anomalous diffusion becomes A ~ — 27ln(f2 cut /f2o)j where f2 cu t is 
a suitable cut off frequency for the Ohmic environment. Thus, the vacuum fluctuations of the environment is felt 
primarily through the anomalous diffusion coefficient that can have a large magnitude. Note that in a macroscopic 
device such as a single Josephson junction biased by a fixed external current 0, IE H EE] one assumes an Ohmic 
environment just to model the junction resistivity. 

Equation (|42Jl is often used to describe the effect of decoherence produced by the diffusion coefficient to study the 
emergence of classical behavior in quantum systems; this is a topic of recent interest; see Ref. J42j for a review. Of 
particular relevance to our problem is the study of decoherence in quenched phase transitions |6l|. and the effect of 
decoherence in quantum tunneling in quantum chaotic systems SI S3 , or in a double- well potential S3 ■ 

The reduced Wigner function W(x,p) describes the quantum state of the open quantum system, and given a 
dynamical variable A(x,p) associated to the system its expectation value in that quantum state is defined by, 



(A(x,p)) = J dxdpA{x,p) W(x,p). (43) 

Then one can easily prove from Eq. I|42(l that defining, 

N = J dxdp W(x,p), {E) = J dxdp (jLj + U(x)jw(x,p), (44) 

we have N = and (E) — ~^((p 2 /M) — Na 2 ). Note that the value of a 2 ~ \Kl for the zero temperature case is 
reasonable since for the virial theorem the average kinetic energy is half the energy, (p 2 /M) = Na 2 , and the averaged 
energy is conserved (E) = 0. 

Further insight in the effect of the different terms of Eq. i|42[) can be obtained from the so called linear entropy 
S = 1 — Tr p 2 , where p 2 is the reduced density matrix of the system. This entropy is also a measure of decoherence 
[62| since for a pure state p 2 — p r and S = 0, whereas for a mixed state Tr p 2 < 1 and S > 0. In position 
representation p r = p r (x,y) and Tr p 2 = J dxdy p r {x,y) p r (y,x). In terms of the Wigner function, see Eq. J27J, 
p r {x,y) = J dp exp[— ipix — y)/h)W[(x + y)/2,p\. If we now call X = (x + y) /2 and u — x — y, then dxdy = dXdu 
and we can write 

Tr p 2 r = 2irh J dXdp W 2 (X,p) . (45) 

We may now compute dS/dt using Eq. (|42|l . The first three terms are total derivatives and do not contribute to 
dS/dt. The dissipation term contributes as 

dSdis 



= -27rn 7 J dXdp W 2 (X,p) , (46) 

which is negative and may be understood as reduction of entropy by heat transfer to the environment. The normal 
diffusion term contributes as 

d l^L = A^M„ 2 I AX A* f^C\ 



dt 



dt -^M^JdXdp^—) , (47) 

which is positive and always increase the linear entropy. The anomalous diffusion term, on the other hand, contributes 
as 

dS ^-^AfdX d p^, (48) 



dt J r dp dX 

which has no defined sign. We may infer from here that normal diffusion will always induce decoherence 
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A. Energy representation of the reduced Wigner function 

Let us now use the base of functions in phase space W Ei e 2 i, x iP)i introduced in Eq. (|32|l . to represent the reduced 
Wigner function W (x, p, t) as in Eq. (|35[) . The previous N and (E) have very simple expressions in this representation: 



N 



J dE C EE (t) , (E) = J dEE C EE (t) . 



(49) 



To check the last equation we note that J dxdp [(p 2 /2M) + U (x)] W Ei e 2 { x iP) = E\5 (E\ — E 2 ) , which can be easily 
proved by explicit substitution of the definition of W ElE2 , and trading powers of p by derivatives with respect to y 
into expressions and partial integrations. 

The quantum transport equation i|42|) in the energy representation becomes, 



d -i 

0j_C ElE2 (t) = — (E 1 - E 2 ) C ElE2 (t) + I dE[dE' 2 Q ElE2 , E > lE > 2 C E ' iE ! 2 (t) , 
where, after one integration by parts, 

dxdp ( d 



(50) 



(51) 



which has the contributions from the dissipative, normal diffusion, and anomalous diffusion parts, respectively, as 

(52) 

From Eq. I].'j2[l it is easy to see that these coefficients can all be written in terms of the following matrix elements: 



!E 1 E 2 ,E' 1 E' 2 — ^ ElE2 ^ E ' lE ' 2 1" ^E 1 E 2 ,E' 1 E^ "T" E 1 E 2 ,E' 1 E' 2 



X 



E X E 2 



Pe 1 e 2 

( XP ) El E 2 



X 



-Ei E 2 



dx X IpE! (x) lpE 2 (x) 



j J dx ip El (x) ^^e 2 (x) , 
Y [ dx xip El (x) -^,i ) E 2 (a?) , 



dx x 2 (x) ip E2 (x) 



Explicitly, we have that 

n (D) ^jl 

^E 1 E 2 ,E' 1 E 2 — 2h 



(XP) EiE , 6 (E 2 - E' 2 ) - P ElE[ X E2E , - X ElE[ P E2E , + (XP) E2E , S {Eh - E[) 



(53) 
(54) 
(55) 
(56) 

(57) 



Q 



(N) 



7 Ma 



2 r 



2X Ei e' 1 Xe 2 e 2 - X EiE ,^5(E 2 - E' 2 ) - X E ^ 6 {Ex - E[) 



(58) 



Q 



(A) 

E 1 E 2 ,E' 1 E! 2 



A 

Ti 



(XP) EiE , S (E 2 - E' 2 ) - P ElE ,X E2K + X ElE ,P E2E , - (XP) E2E , S (Ex - E[) 



(59) 



Thus, in terms of the Wigner function elements C ElE2 the dynamics of the quantum transport equation (|50|l is simple. 
Note that the coefficients and Q( N > preserve parity, while and a coefficient (which corresponds to the 
pure quantum first term in the equation) change parity; in the sense that the symmetric and antisymmetric parts of 
C ElE2 are independently preserved or interchanged when contracted to these terms. 

Equation ijSTJjl resembles a similar equation when a Floquet basis of states are used 13 EE Ell > which are very 
useful when the Hamiltonian of the system is periodic in time. The Floquet basis is discrete in such a case and a 
numerical evaluation of the corresponding matrix elements (|53|I - H56|I can be performed; see for instance |43l l44j| for a 
recent application. It is remarkable that in our case approximated analytic expressions for these matrix elements can 
be found. 
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B. Some properties of the matrix elements 

The matrix elements (|53|) - (|56|l have a clear physical interpretation and several relations can be derived among 
them. Note that X Ei e 2 is the matrix element of the position operator X in the energy representation. Since 
XipE ( x ) = x ^e ( x ) ■ we must have J dE\ Xeei^Ei ( x ) = x i>E ( x ) ■ 

On the other hand, Pe ± e 2 is the matrix element for the momentum operator. The canonical commutation relation 
[P, X] = —iH, implies [H, X] = (—ih/M)P, and taking matrix elements on both sides we have 

( Ei -E 2 )X Ei e 2 = ~^P Ei e 2 . (60) 
Also, X EE is the matrix element of X 2 , therefore 

Xe 1 e 2 — J dE X Ei eX E e 2 - (61) 

On the other hand, (XP) E E is the matrix element of XP, consequently \(XP) E E ] * = — (XP) E E cor- 
responds to PX, and (XP)^ + {XP) EaEi = [X,P] EiE2 = ihS (Ei - E 2 ). Also\xP) EiE2 - {XP^ Bl = 
(iM/K) {Ei — E 2 ) X EiE2 , where the commutator [H, X 2 ] has been used in the last step, therefore 

(XP) ElE2 = ^( E ^ E2 ) X k E2 + j6 (Ei ~ E 2 ) . (62) 

We have, also, that (XP) EiE2 = J dE Xe 1 ePee 2 - One may check, for consistency, that these relations imply N = 0. 
In Appendix B a test of the quantum transport equation in the energy representation (and of the above matrix clement 
properties) is given by checking that a stationary solution with a thermal spectrum is, indeed, a solution in the high 
temperature limit. 

C. Computing the matrix elements 

The matrix elements contain singular parts coming from the integrals over the unbound region beyond x s . These 
singular parts are easy to compute, since far enough the wave functions assume the simple form ©. When performing 
the calculation of the singular parts of the matrix elements we will use that when x — > 00, we have the identities 

sin(px/h) cos(px/h) 

>o(p), ► 0, (63) 

np p 

which can be easily checked by taking the Fourier transforms of these functions with respect to p. 

The computation of the singular parts of the matrix elements I|53|I - (I56|I may be reduced to the evaluation of three 
basic integrals. These integrals are 

Aa,s (PiiPz) = J dx sm[(pi^fp 2 ) x /h + SiTh}, (64) 

and 

B(pi,p 2 )~ J dx sin (pix/h + 5i) sin (p 2 x/h + S 2 ) , (65) 

where, for simplicity, we have written pi = p^ (Ei) and Si = S (Ei) (i = 1, 2). The matrix element Xe 1 e 2 is 

2M f 

Xe x e 2 ~ 7 , dx x sin (pix/H + 5i) sin (p 2 x/h + S 2 ) 

nTTy/prpl J 

f)A f)A /r)S, rM„\ 

(66) 



A I 



9A__dA__ fdSi 862 
dpi dp 2 \dpi dp 2 



B 



where A = (A s - A A )/2 and A = {A s + A A )/2. The matrix element X% Ei is 

X EE ~ - — ■ / dx x 2 sin (pix/h + Si) sin (p 2 x/h . + S 2 ) 

- - h-K^/piPi J 

2M 



K^JP\P2 



dC _ f dSi 
dpi [dpx 



D 



(67) 
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where, it is easy to show that C = (dB/dpi) — (dSi/dpi)A, and that D = —(dA/dpi) — (d8\/ dp\)B . The matrix 
element Pe^e 2 is 

P Ei e 2 ~ z ; (Pi +P2 A 68 

which according to the relations among matrix elements derived in the previous subsection is related to Xe 1 e 2 by 
Eq. H60|) . The remaining matrix element (XP)e 1 e 2 , on the other hand, can be computed from the element X^ iE2 
according to Eq. 

1. The integrals A(pi,p2) and B(p\,p2) 

Thus, we are finally left with the computation of the integrals (|64f) and (|65|l . The integral B (pi,£>2) of Eq. (|(j5|) is 
dominated by its upper limit x 



1 f s 

B{pi,P2) ~ 0/ dx C0S ^ Pl ~P2)x/h + 5 1 - 5 2 ] 



sin [(jpx - p 2 ) x/h + Si- 5 2 ] 



2 (pi - J*) 



7T^ 

-y^CPi-Pa), (69) 

The integrals Aa,s (piiPz) are more subtle. The integral As is clearly regular on the diagonal. Since we are interested 
mostly on the singular behavior of the matrix elements, we can approximate ^5 ~ 0. On the other hand Aa is 
exactly zero on the diagonal. Close to the diagonal, the integral is dominated by the region where the argument of 
the trigonometric function is small, and thereby the integrand is non oscillatory. Estimating the upper limit of this 
region as x ~ h {p 1 — p 2 ) 1 , we get 



A A - h~ L (pi - pa) x 4 + (Si -5 2 )x = hPV + ... , (70) 

\Pi -Vi) 

where the dots stand for regular terms. Actually, this argument would allow us to introduce an undetermined 
coefficient in front of the principal value PV, but in the next section we show that h is the correct coefficient, as 
follows from the canonical commutation relations. 

Thus, we are now in the position to give the explicit expressions for the singular parts of the matrix elements and 
write, finally, the quantum transport equation in its explicit form. This is done in detail in the next section. 

V. THE QUANTUM TRANSPORT EQUATION 

In this Section we explicitly compute the quantum transport equation (|42(l satisfied by the reduced Wigner function 
in the energy representation. 

A. Matrix elements 

First, we need to compute the matrix elements described in section TlV CI We begin with the matrix element Xe ± e 2 
which according to l|66(l and (|69|I - H70(I can be written as: 



Mh 

■^E X E 2 



1 ".pvC-L-V*.^-^ 



7r dpi \pi — p 2 ) dpi 



VP1P2 

We go next to the matrix element PexE 2 , which from l|68f) and H7U|I can be written as, 



(71) 



P Ei e 2 = -^=1-(pi+P2)PV^— . (72) 

y/PlP2 27T pi - p 2 
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These two operators X and P are connected through Eq. i|60|) . It is easy to check that the two previous results satisfy 
this relation. Just notice that from Eq. (|A17|I we can write E\ — E 2 = (p\ — p 2 )/2M which together with Eq. (|71|) 
for Xe 1 e 2 l ea< i to —ih/M times the right hand side of Eq. I|72|l. that is 

{Ei - E 2 )X Ei e 2 = --^Pe 1 e 2 - 
Another check of the previous results is the consistency with the canonical commutation relations 

J dE {P Ei eX E e 2 - X EiE Pee 2 ) = ~ihS (E 1 ~ E 2 ) . (73) 

This check requires a little more work. First it is convenient to change to momentum variables and write, S (E\ — E 2 ) = 
(M I y/pip 2 )5 (pi — p 2 ) . Then one needs to compute the integral 



i = h 



dpPV 



1 



Pl-p 



PV 



1 



P-P2 



= -hn 2 6 (pi - p 2 ) , 



(74) 



The evaluation of this integral is easily performed using the following representation of the principal value 

PV {1) = / ^^H^gnKD, 

which is easily proved by taking the Fourier transform of PV(l/p). With the result of Eq. I|74l) it is straightforward to 
check that the commutation relation (|73|l is an identity within our approximation. This consistency check is important 
because it can be used to fix to h the coefficient in front of the principal value of Aa in the argument leading to Eq. 

G3- 

We can now move to the matrix elements for X 2 . Having an expression for Xe 1 e 2 i n Eq. I|71|l it is best to compute 
X\ E directly from the relation l|61|) which leads to 



Y 2 

EiE 2 



2 r 



Mh 
VPTP2 

'dSi 
.dpi 



9 2 s , \ f (dSx 

(Pl - P2 



dpidp; 



dS 2 



7r \ dpi dp 2 J dp 2 







-PV 



Pl ~P2 



8 (pi - Pi) + 



(75) 



where we have used the result ()74|l and performed the E integration or, more precisely, the p integration. 

The matrix element (XP)e 1 e 2 = J AEXe 1 ePee 2 can be analogously obtained from the expressions l|71|l and 17211 . 
The result is 



(XP) 



EiE 2 



iMK 



d 1 dS ( 1 

2 P2^— S(pi -Pi) + -7T 1 (Pl +P2)PV I 

dpi 7r dpi \pi - p 2 



(76) 



A further consistency check of these expressions comes from the property H(j2|) . which is satisfied within our approxi- 
mation. 



B. The quantum transport equation and time scales 

Finally, we can write the quantum transport equation (|50|l . in a more explicit form. The coefficient Q is given 
by (JH21, with the values of the dissipation and diffusion parts given by fB!^ and fBl^l . which can be directly 

computed using the matrix elements obtained in the previous subsection. It is convenient to introduce new Wigner 
function coefficients, 

M 

C Ei e 2 (t) = C PlP2 (t) , (77) 

VP1P2 

and the result is the rather cumbersome expression (|C1|) given in Appendix |0 As explained there we can get a local 
approximation of the quantum transport equation <|C1|) : 

dC(P,p,t) f-iPp <9 „, 2 d 2 3\ , n , 

- 1 +l^P + jMa 2 —+iAp—)C(P,p,t) 



dt \ Mh ' dP ' dP 2 r dP 

.2 fd(Si-s 2 ) , d(Si + s 2 y 2 
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where the average and difference momentum variables P and p, defined in Eq. (|C3(1 . have been used and where 
C(P,p,t) = C piP2 (t). From the quantum transport equation (|C1I) . or its local version Eq. (|78|) . it is easy to discuss 
the different time scales of the problem. The first, of course, involves the dissipation term which includes the dissipation 
coefficient 7, it defines a time scale ~ 7 _1 which is the relaxation time. 

But before we go on with the interpretation of the different terms, it is important to recall the meaning of the 
Wigncr function coefficients C PlP2 , or Ce 1 e 2 - First, we note that these coefficients are directly related to the co- 
efficients Ce of the energy eigenfunctions which make the tunneling state from the false vacuum in the isolated 
system, i. e. when there is no interaction to the environment. Thus, the coefficients Ce 1 e 2 describe the quantum 
correlations between the energy eigenfunctions that make the tunneling system. These coefficients are initially sep- 
arable Ce 1 e 2 (®) — Cbi(0)C|, 2 (0). In the isolated closed system its time evolution, as given by Eq. (|3l)|) . is simply 
Ce^e-i (t) = Ce 1 e 2 (0) exp[—i(Ei — E 2 )t/h] , which means that these correlations keep their amplitude in its dynamical 
evolution. 

This is very different in the open quantum system as a consequence of the negative local normal diffusion term in 
Eq. (|C1J) which depends on the phase shift derivatives, i. e. the last term of Eq. Ij78(l . or equivalently the term (|C2|I 
when it is written in the p\ and pi variables. This negative defined term has no effect for the diagonal coefficients, 
when Ei = Ei, but its effect is very important for the off diagonal coefficients. In fact, the amplitude of the off 
diagonal coefficients exponentially decays in time. The time scale can be estimated by taking the derivatives of the 
phase shifts 5i (i = 1,2) near the false vacuum energy Eq, which is where the energy wave functions pile up. Using 
Eq. (|T2l it is easy to see from expression <|C2|) that this time scale is of the order of 

(79) 

where r# is the relaxation time, Xb = h/ (2ayM) is a characteristic de Broglie wavelength (in the high temperature 
case when a 2 = ksT it corresponds to the thermal de Broglie wavelength), and Id ~ a 2 h\/Eo + Uoo/(eyM) is a 
characteristic length of the problem with a a dimensionless parameter that measures the scale of the energy differences 
of the off diagonal coefficient, Ei — E2 ~ as; so it is of order 1 when the energy differences are of order e. Thus, the 
last term of equation (|78l) destroys the quantum correlations of the energy eigenfunctions. The time scale td may be 
considered as a decoherence time [3^ , and thus the effect on tunneling of this term may be associated to the effect of 
decoherence. 

Another time scale in the problem is, of course, the tunneling time which according to (|26|l and © is given by 
Ttunn ~ h/e. Its relation to Td is given by Td ~ itunn/ (o^-D), where the dimensionless parameter D is 

It seems clear that when te "C T tunn the coefficients Ce 1 e 2 become diagonal very fast and the local approximation 
to the transport Eq. (|78|l is a useful approximation. 

VI. TUNNELING IN THE OPEN QUANTUM SYSTEM 

We can now compute the tunneling rate from the false vacuum for our open quantum system. Thus, let us assume 
that our particle at t = is trapped into the well of the potential in the false ground state with the energy Eq, 
i. e. the ground state of the auxiliary potential U aux introduced in Section III CI We know from that section that the 
wave function tpo of this state can be expressed in terms of the eigenfunctions ^ by Eq. I|13|) with the coefficients 
Ce given by Eq. 

Under the assumption that the decoherence time is much shorter than the tunneling time the quantum transport 
equation (|C1|) may be simplified. After a typical decoherence time the Wigner function coefficients Ce ± e 2 become 
diagonal. We will begin our leading order approximation assuming that these coefficients are diagonal, which means 
that decoherence is almost instantaneous, and then we will correct it in a perturbative way. 

A. The Kramers limit 

We are interested in the regime where the decoherence time is much shorter than the tunneling time, or D > 1 i. e. 
when the decoherence term (|C2|) suppresses the non diagonal Wigner function coefficients C PlP2 (t) and get a totally 
decohered Wigner function. In this case using the average and difference momentum variables P and p, see Eqs. <|C3ll . 
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we may approximate these coefficients as C PlP2 (t) = C(P,p,t) ~ f(P,t)S(p). The quantum transport equation I|C1|I . 
or its locai approximation Eq. J7SJ), reduces then to 

< 81 > 

which is the Kramers [17| thermal activation equation for f(P,t). This is a continuity equation for a distribution 
/ with a flux $ = —7 [P + Ma 2 J=p) f . The computation of the tunneling amplitude has thus been reduced to the 
computation of the escape probability of a particle confined to a potential U(x) defined in Eq. subject to a 
damping force 7P and white noise with amplitude jMa 2 ; which corresponds to a temperature T — a 2 /ks- The 
boundary conditions of Kramers problem are the usual ones. We assume that the particle is initially trapped at the 
potential well. Let us define Po = y / 2MU 00 , we may take Uoo = in this section since the asymptotic value of the 
potential outside the barrier has no relevance for the classical activation problem. We will ask that there is no flux 
entering into the well so that $(0) = and df/dP(0) — 0. Moreover, at the separatrix when P s = y / 2Me s , where e s 
is the potential barrier, f(P s ) = 0. This means that the number of particles above the separatrix is negligible. 

To find the activation rate we seek normal modes f = f (P) e~ rt and assume that r is very small. If r = we have 
the stationary solution f — exp(— P 2 /2Ma 2 ) which satisfies the boundary condition at P = but not at P s . We seek 
a second stationary solution, for r = 0, of the form Fo = fof s then Eq. (|H l|t reduces to 

lMa2 ^p{ f °Jp fi ) = °' (82) 

which leads to df s /dP — K/ fo, where the value of the constant K is irrelevant, and may be chosen as K = — 1. The 

resulting solution Fq(P) — fo J p s dQf ~ 1 (Q) satisfies the boundary condition at P s but not at P = 0. Thus, we now 
have two stationary solutions of Eq. 181|) . fo and Fo, and we may use the variation of constants method to find a 
normal mode solution for r > 0. Let us write 



f = a(P)f +P(P)F , 



(83) 



with the supplementary condition a' fo + fl'Fo = 0. The function a{P) and f3(P) are determined by substitu- 
tion into Eq. I|81[) which leads to the first order differential equations a' = — (r/jMa 2 ) (aFo + (3Fq / fo) and 
P' = (r/jMa 2 ) {afo + (3F ); the boundary conditions are 0(0) = a(P s ) = 0. For r very small, we may just ap- 
proximate a = a (0), and (3 = in the right hand sides of those equations and we get 



a(P) = a (0) 1 



dQ Fo (0) 



(84) 



which leads to the lowest eigenvalue 



7MB0 



Jo 8 dQ Fo (Q) 



(85) 



This is the escape rate of the particle. We may give this rate in terms of the effective escape temperature defined in 
Eq. by equating r to (l/2r) exp(— e s /fcsT esc ), that is 



T 



a 

k B 



■In 



27 fne 



(86) 



where we have used that the dynamical time r, defined in Eq. 
of the order of T esc ~ a 2 /ks < ^Mlo/ks, at zero temperature. 



Oo V 

JH1, is r ~ 7r/57 . Thus, we have a escape temperature 



B. Effect of anomalous diffusion 



We may now estimate the effect of anomalous diffusion into the tunneling rate. We have seen in section Hvl that 
the effect of anomalous diffusion to the linear entropy of the reduced density matrix of the open quantum system is 
undefined, unlike the effect of normal diffusion. It may increase or decrease the entropy, and hence the decoherence, 
depending on the product of gradients of the Wigner function on phase space; see Eq. i[48|l. Let us go back to the 
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quantum transport equation (|50|l . The term Q^ N ^ is dominated by the decoherence term l|C2fl : for this reason we 
may distinguish a fast and a slow dynamics. The fast dynamics corresponds to the decay of the non diagonal Wigner 
function coefficients. The slow dynamics is the diffusion of the Wigner function coefficients along the diagonal which 
will be considered in the next section and may be described by a Fokker-Planck equation for a classical distribution 
function. 

To evaluate the effect of Q^ A \ let us consider Wigner function coefficients Ce 1 e 2 of the form 

C = C diag + C odd , (87) 

where C dlag is diagonal and C odd is antisymmetric. Here we do not include symmetric off diagonal coefficients, such 
as C even , because we assume that decoherence takes them to zero, even though other terms generate it. According 
to the parity properties of the Q^ D \ and defined in Eqs. (|E7jl . lt55|) and we may write Eq. (|l)0|) as a 

set of coupled equations: 

Sfjdiag = (q{D) + q(N)\ (jdiag + q(A) (jodd /ggx 

at V / 

^- t C odd = (Q (D) + Q (Ar) ) C odd + Q^C diag . 

Note that the first term of Eq. (|50|l does not contribute to this set of equations: Q^' does not contribute to the 
first equation because it vanishes on the diagonal, nor to the second equation because it destroys diagonal Wigner 
functions. Note that the symmetric off diagonal coefficients C even do not couple to the diagonal coefficients C d%aa 
through terms that contain a 5{p), thus these terms are second order with respect to the diagonal coefficients as 
analyzed in section TVI Al Because the leading process is decoherence, the second equation may be approximated by 
d_Qodd = =±C odd + Q^C dla9 where td is a typical decoherence time scale estimated in section IVBl see Eq. (JJ5). 
Note that this is analogous to the relaxation time approximation for the Boltzmann equation. If td is short enough, 
C odd simply trails C dlag and C odd ~ tdQ^ C dia9 and we obtain an autonomous equation for C dlag 



8 

radian 

d~t 



(jdiag _ (gg) 



To compute (Q^) 2 we go back to Eq. Q59[l. we permute the operators X and P so that we get an expression 
in terms of P and P 2 only. The linear terms in P vanish near the diagonal and then replace the P with a typical 
momentum scale, such as v Ma 2 . Finally we get 



( Q(A) ) 2 Ei e 2 ,e IE > ~ {xkE'MK + 6 ElEl X 2 2E , - 2X ElE ,X E2E ,) , (90) 



which is the same operator as . Thus, equation (|89|l leads to the first of the previous set of coupled equations 
when the anomalous term is neglected, but with a modified normal diffusion coefficient. The effect of anomalous 
diffusion is then to lower the normal diffusion coefficient a 2 to c^, where 

This translates into a lowering of the effective escape temperature of Eq. (|86p. namely T esc ~ cr 2 ff/ks- Thus 
the escape temperature is always lower than a 2 /ks but it tends to it when the decoherence time td goes to zero. 
Consequently, the overall effect of anomalous diffusion is to inhibit tunneling when strong decoherence is assumed. 



VII. CONCLUSIONS 



To conclude, let us briefly summarize our results. Under the assumption of strong decoherence and using a real 
time formalism we have estimated the tunneling rate for an open quantum system representing a quantum particle, 
trapped in a local minimum of a cubic potential, coupled to an environment. The real time formalism is based on 
the master equation for the reduced Wigner function that describes the open quantum system. Our computational 
method involves the introduction of an energy representation of the reduced Wigner function which is based on the 
energy eigenfunctions of the isolated system. The master equation in this representation, Eq. ||SUJ), is an equation for 
some Wigner function coefficients that describe the quantum correlations between eigenfunctions of different energies. 
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In spite of its apparent straightforwardness, the original master equation i|42[l contains a term with three derivatives 
of the Wigner function, which makes it quite hard to handle cither analytically or numerically. From this point of 
view the reduction of Eq. (|42|l to Eq. H5()(l . which is then approximated by Eq. (|78|l . is already a definite step forward. 
Also, the basic process of decoherence is hard to discern from Eq. (|42[1 . while it is clearly depicted in Eq. I|78l) . A 
striking feature of the decay process is that in last analysis it is robust against the details of the potential barrier. 
For example, the Kramers activation rate is dominated by a single feature of the barrier (its height), and similarly 
the quantum tunneling rate for the closed system. This simplicity is lost in Eq. I|42|) but regained in Eq. (|78|l . We 
see that the details of the potential (such as the height of the barrier) enter the problem only through the boundary 
conditions and the phase shifts in the energy eigenfunctions. In particular, the strength of decoherence is determined 
by the phase shifts near the false vacuum energy, where the simple approximation of Eq. H12|) applies. 

In our problem, where the particle is initially trapped in the false vacuum, the master equation involves a term, ljC2fl . 
that destroys the quantum correlations of the eigenfunctions and is, thus, responsible for decoherence. The strength 
of this term is characterized by the dimensionless parameter D, defined in Eq. I|80|l . which is directly proportional 
to the energy difference between the false and true vacua. Under the assumption of strong decoherence the pure 
quantum channel to tunneling is partially suppressed; since decoherence destroys the fine tuning among the energy 
eigenfunctions that makes tunneling possible in the isolated closed system. Tunneling then follows an activation- like 
channel due to the zero point fluctuations of the quantum environment. This is similar to the result recently found 
in Ref. |45| for a double well potential, which uses a large scale numerical simulation to solve the master equation. 

Thus the picture we have is the following. For the isolated closed system tunneling from the ground state goes 
through the usual quantum channel which the WKB approach or the instanton approach reproduce, and that we 
can equally reproduce using the energy representation of the master equation; see section 1111 Bl When the system 
is coupled to an environment the general effect is felt as dissipation and diffusion, the latter coming as normal and 
anomalous diffusion. These terms will produce essentially two effects. On the one hand they will produce decoherence 
to the system which will tend to suppress tunneling as the system becomes more classical. On the other hand diffusion 
will also introduce noise into the system which will induce tunneling by a mechanism similar to thermal activation. 
In general all these mechanism act simultaneously and their effects cannot be disentangled. 

When decoherence is very fast quantum tunneling is strongly suppressed and activation dominates, to leading 
order. In section IIVI we have seen how the diffusion term, the normal diffusion and anomalous diffusion affect the 
linear entropy of the reduced matrix density of the open quantum system. We see clearly that normal diffusion tends to 
always increase the entropy and thus to induce classicalization and decoherence. This is in fact what the decoherence 
term (|C2|) in the quantum transport equation does. The fact that activation from normal diffusion becomes the 
dominant decay mechanism when decoherence is fast enough validates the analysis of vacuum decay in cosmology and 
quantum field theory given in Refs. [25l l33l l34|. On the other hand, these references show how the analysis given 
here may be improved by a more realistic description of the bath and the system-environment interaction. 

We should emphasize that although the master equation in the energy representation, Eq. l|5Uf) . is much simpler 
than in the standard phase space representation, Eq. l|42ll . it may be difficult to go beyond the present results by 
analytic means. It may still be possible to treat perturbatively the symmetric but non diagonal terms in the Wigner 
function coefficients. These terms were neglected in section lYl Bl where we considered the effect of anomalous diffusion 
assuming coefficients of the form (I87f) . The quantum coefficient which induce pure quantum tunneling in the 

isolated system, will have an imprint in those terms. 
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APPENDIX A: WKB SOLUTION 

In this Appendix we solve the WKB problem posed in section Hi Bl The starting point are Eqs. ©, Q and © 
with the cubic potential of Eq. @, we have to match the WKB solutions in the different regions across the potential 
function. 
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1. Matching from forbidden to allowed regions 

Let Xo be a classical turning point U (xo) — E, and let U' (xq) < 0. Then to the left of xo we have a forbidden 
region, the two corresponding independent WKB solutions of the Schrodinger equation J3J are 

e ±S(x ,x)/h 

F±(x ,x) — — (Al) 
V2p(a;) /ft 

whereas to the right of xq the two independent solutions are 

G±(x,x ) = —===, (A2) 
^2p (x) jh 

and we wish to find the corresponding matching conditions. For x — > Xq~, we can Taylor expand the potential around 
xo and write p (x) = k (xo — x) 1 ^ 2 and S'(xo,x) = |k(xo — x) 3 ^ 2 , where we have introduced K = y/2M \ U' (xo)|. 
Similarly for x — > Xq, we have p (x) = k (x — xo) 1 ^ 2 and S (x, xo) = |fc (x — xo) 3 ^ 2 . 

If we write x — Xo = e™ (x — x) then iS (x, x ) = S (x , x) and it would seem that simple analytical continuation 
yields G+ (x,xq) — > e~ l7r ^F + (xo,x). However, this is impossible, recall that if we define the flux J = —i(ip*d x ip — 
ipd x ip*) then the Schrodinger equation implies flux conservation d x J = 0. Now G+ (x, xq) has J — 1 and therefore it 
cannot turn into F + (xo,x) , which is real, and has J = 0. Thus, we try instead 

G+ (x, x ) - e- <7r / 4 F + (x , x) + /3F_ (x , x) , (A3) 
and imposing flux conservation we obtain /? = (1/2) exp(i7r/4). We therefore find the matching conditions 

e T-/4 F+ (a . 0) x) + l e ±ftr/4 F _ ( XOj a.) ^ G± (!C) Xq) t ( A4 ) 



from were we finally obtain, using Eq. (|A2I 



F + (x ,x)^ = — cos [ ^5(x,x ) + 7 ) , (A5) 



V 2 P W A 



a ~V»~u, 4 



and 



F_ (x ,x) -> 7==== sin ( (x,x ) + j ) ■ (A6) 
v/2p (x) /h \h 4J 

2. Matching from allowed to forbidden regions 

Now consider the case when U' (xo) > 0. To the left of xo, we have an allowed region and the solutions are oscillatory 
G± (xo, x), to the right of the turning point we have a forbidden region and the solutions are a linear combination of 
(|A1I) . By exactly the same procedure of the previous section, after imposing flux conservation across Xo we obtain: 

G T (x , x) -» e***l*F + (x, xo) + \e^^F^ (x, x ) . (A7) 
Note from these equations that the solution that matches a decreasing exponential is 

' :sin (-S(x ,x) + ^ ] -F- (x,x ). (A8) 



^2p (x) /h V s V 2 

3. WKB solution for < E < e s 

We can now put all this together to write the energy cigenfunctions for our cubic potential @ for energies in the 
range < E < s s . There are three classical turning points in this case xl < xr < x ou t- To the left of xl we have a 
forbidden zone extending to -co, so we have 

ipE (x) ~ K E F- (x L , x) ; x < xl, (A9) 
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where Ke is a normalization constant to be determined latter. To the right of xl we have from (|A6|I 

fcw ~vs§b s Ks S( " i,+ £>' (A10) 

which after using the definition JSJ can be rewritten in the region xl < x < x R as 

IpE (X) ~ K E ^(S(x R ,x L )/h-^/i) G _ ^ x) + e -i(S(x R ,x L )/fi-7r/4) G+ ^ ^ _ (All) 

This expression is in the form suitable for extension to the forbidden region, that is, to the right of x R . Thus, by 
using l|A7(l we have to the right of x R 



ip B (x) ~ 2K E 
which can be rewritten again as 
ip E (as) ~ 2K E 



cos ( -S (x R ,x L )j F + (x,x R ) + isin (-S {x R ,x L ) ) F_ (a\. <■/,•) 



(A12) 



which is in a form suitable for extension to the right of x ou t' 

2K E 



(A13) 



tp E (x) 



■2 <-, .s ( is (a*, x L ) J e S(x«.«^)/» S in f is (x, x out ) + J 



y/2p{x)/h 

+ i sin Qs (a*, a*)) e -S(*~*.*i«)/» C o S Qfi ( X) x ou( ) + 1 



Note that if we impose the Bohr-Sommerfeld quantization rule 

7T 

S (x R ,x L ) = -(1 + 2ra)7i, 



(A14) 



(A15) 



only the subdominant, exponential decreasing part survives. This would correspond to the case when the far right 
region is forbidden and may be used to define energies for false states trapped into the potential well, in particular 
n = will correspond to the false ground state. 



4. Normalization 



All that remains now is the determination of the normalization constant Ke which can be done from the normal- 
ization of the wave functions. The eigcnfunctions are subject to continuous normalization 



dx tpE! (x) ipE 2 (x) = 5 (Ei - Ei) . 



(A16) 



Since the functions themselves are regular, the singular behavior must come from the upper limit, see for instance 
|47j . For large enough x, we have from Eq. © 



Poc = ^2M{E + U oc ). 



Let us write from Eq. tKl\ . 

S (x, X ou t) = Poo \X — Xgut) + 



dx' 



y/2M {E-U (x')) - y/2M (E + U x ) 



if this integral converges, we may take the upper limit of integration to oo, whereby 

S (X, Xout) = PooX + f (E) , 

where f(E) stands for the second term of (|A18II . Then, for x ^> x ou t, we can write from (1 A 14(1 and l|A19[) 

/2hK E 



IpE (x) 



/Po 



A (E) sin (^f) + B (E) cos (^) 



(A17) 



(A18) 



(A19) 



(A20) 
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where A(E) and B(E) are given by 

A(E) = 2cos(j:S (x R ,x L )^j e s ^ x «V h cos (f(E) + J) 

-^sin (^S(x R ,x L )^ e -s(.«*.» ^sin (/ (£) + J) , 



(A21) 



= 2cosQs , (x fl ,a; L )je' s ( ;l; °" t ' :Efi )/ ?i sin(/(£;) + £ 

+i sin (a*, a*)) e -S(xo.t,x«)/ft cos ( ; (s) + ^ . 
Substituting Eq. l|A"20|) into JXTSJl, the singular terms in the normalization integral are 



(A22) 



K 



dx $ El (X) *Pe 2 (X) ~ h 2 TT— ^ [A 2 + B 2 (Si)] 

Pi 



dp 



dEi 



-l 



5 - Ea) , 



where the delta function comes from the x integration which brings S(p\ ~pi) and where we have defined pi = Poc(Ei) 
(i = 1,2) and changed from momentum to energy variables according to Pidpi — MdEi\ see Eq. (|A17|) . The 
normalization condition reduces to 1 the coefficient of the delta function above 

tfn^k [A 2 (E 1 ) + B 2 (E x )] = 1. (A23) 

This suggests the introduction of the phase 5 E as follows, 



K E A(E) = x /-M-cos5 E , K E B(E) = ^^-smS E . (A24) 



Thus, the eigenfunction at i> x ou t is Eq. ©, that is 



f \ I 2M ■ P 



ipB (x) ~ \ j sin (E^^t + § E 



To work out the constant K E in greater detail we note that form Eqs. I|A21() and (|A22|) we have 

A 2 + B 2 = 4 cos 2 Qfi iz)) e 2 ^(— + I sin 2 Qs (a*, x L )j e~ 2S ^ x ^ h , (A25) 

which is non vanishing as long as E is real. However, if we allow for complex energies, as is typical of unstable states, 
it may be zero provided 

cos 2 (Is (x R , x L ) \ = ^ sin 2 (Is (x R , x L )\ e -^(*»*.*i0/» (A26) 
The left hand side is zero whenever the energy satisfies the Bohr-Sommerfeld condition (|A15|) . 

APPENDIX B: THERMAL SPECTRUM 

In this Appendix we check that the quantum transport equation l|50() admits a stationary solution with a thermal 
spectrum. This can be seen as a test on the restrictions satisfied by the matrix elements H53|) - (|56J) with er 2 = fc^T. 
An unnormalized thermal density matrix in the position representation reads, 



p (x, x 1 ) = J dE e~ &E ^ E (x) ip E (a/) , (Bl) 
where (3 = (fc^T) -1 and its associated Wigner function is 

W (x,p) = J dE e~ 0E W EE (x,p), (B2) 
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which in the energy representation in the base We 1 e 2 of Eq. (|35|l corresponds to the coefficients Ce 1 e 2 
e -f3 E i§ (E 1 — E 2 ). Inserting this into the transport equation we get 



J dE e~ f3E Q El E- 2 ,EE = 0, (B3) 
which after using Eqs. (|52p. (|57|) and 158|) can be written in operator language as 

= -LUxPe-PB-e-PBpX-Pe-PBX + Xe-Wp) 
2MH y ' 

+ _L (x 2 e^ H + e~^X 2 - 2Xe-P H X) . 

At the infinite temperature limit, (3 = 0, this is 

= — - (XP - PX) - 4 (^ 2 ff + # AT 2 - 2XHX) . (B4) 
M n n z 

The first term is the commutator which gives — M _1 , and the second term can be written as — -p- [X, [X,H]], which 
using [H, X] = (h/i)(P/M) is easily seen to cancel the first term. 



APPENDIX C: QUANTUM TRANSPORT EQUATION 

Here we write explicitly the quantum transport equation l|50|) in the energy representation. The coefficient Q in 
Eq. H50JI is given by Q52JI. and the values of the dissipative, normal diffusion, and anomalous diffusion parts of this 
coefficient are given, respectively, by Eqs. (|57|l . (|58|l and (|59|l . These parts can be directly written using the matrix 
elements deduced in Section^ When the coefficients C PlP2 defined in Eq. (|77|) are introduced the transport equation 
becomes, 

*9CpiP2 ^ / 9 9\ ^ t. r 9 ( ^ 9 



( Pl - p 2 )C PlP2 + -yMa [ —~ + —J C PlP2 



8t 2MH^ L w-rij* ■ y dp 2 • g p 2 

7 ( d d \ _ . . / 9 9 

7-i2A a 



47T 2 9pi 

7 + i2A 9 

47T 2 9^2 

2 7 Mcr 2 <9 2 



M d P2 (P2 + Pa) P (P2 - P 2 )P(Pl ~ Pl) C p' lP ' 2 

dp' x dp' 2 (pi + p' x ) P(p 2 - p' 2 )P{ P i - piJCp/pj 



■ g pi g p2 J d Pl d P2 P (Pl - Pl) P (P2 - P 2 ) C p > lP > 2 



2 ( dSi 85 2 , 

7 / 9(5i 9(^2 

47r \9pi dp 2 

iA (d8y dS 2 

2ir \dpi dp 2 

J it \dp 2 dp' 2 dpi J dp 2 



dp [(pi + p) P(pi - p)C PP2 - {p 2 + p) P{p 2 - p)C PlP ] 
dp [(pi + p) P(pi - p)C PP2 + (p 2 + p) P(p 2 - p)C PlP ] 



pip 2 



iw- 2 f dp 1 , fdS 1 6S V dS 2 \ dP( Pl -p' 1 ) 
-fMa I —[— + — -2—J — C P{P2 , (CI) 



where we have used the shorthand notation P(x) = PV(l/x). 

This equation may be considerably simplified by noticing the effect of the negative defined local term which depends 
on the phase shift derivatives: 
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This term has no effect on the diagonal Wigner function coefficients, when p\ = P2, but it exponentially reduces the 
off diagonal coefficients C piP2 on a time scale of the decoherence time, t d , as discussed in Section |V] This suggests 
the following local approximation to the transport equation (|Cip whenever tjj <§; T tun n- 

To derive the local approximation it is best to introduce new average and difference momentum variables, 

P =^(Pl+P2), P=Pl~P2,, (C3) 

then when substituting into Eq. IjCljl we have two typical terms involving the p\ and p' 2 integrations: 

J dUdu (2P + U)P(U-u)P(U + u)C(P + U,p+2u), J dUdu (p + u) P(U - u)P(U + u)C (P + U,p+ 2u) , 

where we have used C(P,p) = C P1P2 , and U = P' — P and u — hip' — p). Now we make the hypothesis that 
the P dependence is softer than the p dependence, so within these integrals we can approximate C (P + U,p + 2u) ~ 
C (P,p + 2u). Then using Eq. dJ> we can integrate / dU P(U-u)P(U+u) = ^-S (u) and / dU UP(U-u)P(U+u) = 
0. In this way all the terms in Eq. IClfl which do not depend on the phase shifts Si and 82 become local and 
considerably simplified. For the phase shift terms we may notice that far from the resonance all terms are negligible 
and close to the resonance the local term (IC2(I is clearly dominant so it makes sense to keep only this term. Finally, 
Eq. (|78|l follows as the local approximation of the quantum transport equation IjClll . 



APPENDIX D: TUNNELING RATES FOR THE CLOSED SYSTEM 



In this Appendix we review the calculation of the quantum mechanical tunneling rate for the closed system, that 
is, ignoring the interaction with the environment. The quantum tunneling rate as given by the instanton calculation 

is 



r (inst) _ 0>q _A mil 
1 closed — 2 T ' 



where r = ir/flo, A = Ss/fi = 18e s /(5eo) ~ 12.376, and the prefactor a q — (^OttAo) 1 / 2 ~ 68.306; here we use 
as an example the values e s /fc_B ~ 589.74 mK, and eo/^s ~ 171.55 mK, from a tunneling experiment for a single 
Josephson junction between two superconducting electrodes biased by an external current reported in Ref. 2]. With 
these values, the escape temperature defined in Eq. H41fl is 

Tt St) = ^ 72.345 mK. (D2) 

3.6 - {e /s s ) \na q 

It is interesting to check that this result agrees with the result we obtain when the dissipation is zero. We can use 
our WKB result as obtained in Sec. |nj see Eq. J2EJ), to write 

V {WKB) _ J_ A (m] 
1 closed ~ 2t ' 

where A = (2/h)So{x out , xr), with So defined in Eqs. 0)-©, where the potential U(x) is given by Eq. (0). 

For a cubic potential, the relationship among the energy E, the frequency Q and the action S (xr, Xl) is best given 
in parametric form, 

E = 2e s ((k), n = fio/(*0> S(x R ,x L ) = ^-F(k), (D4) 

"0 



with < k < 1, and 



1 f (1 + k 2 ) (1 + fc 2 ) 1 
C(k )= \2 + 3± j^-A >\, D5) 

8\ [Q(k)] 1/2 [Q(k)f 2 



f(k) = ^[4Q(k)f i K[k^ 1 



5/4 



{a (fc) E [k 2 ] - (1 - k 2 ) b (k) K [k 2 ] } , (D6) 
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where E [fc 2 ] and K [fc 2 ] are the complete elliptic integrals, and we have introduced the functions Q (fc) = 

(1/4) (1 + 14fc 2 + fc 4 ) , a (fc) = (16/15) (2 - fc 2 ) 2 - (1/5) (l - fc 2 ) (21 - 5fc 2 ) and b (fc) = (8/15) (2 - fc 2 ) - (l - fc 2 ) . 

The Bohr-Sommerfeld condition Eq. (|A15|I for the ground state (n — 0), corresponds to the parameter kcs such 
that F (kcs) = Keo/e s which implies that kcs ~ 0.1152. This corresponds to £ (kcs) ~ 0.1423 and / (kcs) ~ 0.9550, 
while the harmonic approximation for the potential yields 0.1454 and 1, respectively. 

To compute the barrier penetrability, A = (2/K)So(x ut, %r), we observe that S (x ut,xn) at energy E is equal to 
S(xr,xl) at energy E re j — e s — E. The exchange of E by E re f is equivalent to the exchange of fc by fc re /, where 
C (fc re /) = 1/2 - C (k) ■ For k G s we obtain k ref ~ 0.2433 and F (k ref ) ~ 2.4073. Therefore 

A = —F (k ref ) ~ 8.459. (D7) 

This is to be compared against the instanton exponent Ao — lna 9 ~ 8.152. In terms of the escape temperature, the 
WKB approximation yields 

T - m = wu s ( £o !% (o 7?n ~ 70 - 869 mK ' (° 8 ) 

F (k r ef) - (eo/£s) m (ilcs/^o) 

which is in good agreement with the instanton result; here we have again used the previous numerical results for e s and 
Eq. This agreement, of course, should not be surprising since for a closed system our method reduces to the standard 
WKB calculation. The purpose of this exercise is just to check the consistency of our calculation and to illustrate 
how the instanton and WKB methods compare. That the difference between exp(Ao) and exp(A) is accounted for by 
the prefactor a q of Eq. IjDlfl can be seen analytically by a perturbative calculation. 
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